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We present the first observation in hadronic collisions of the electroweak production of vector 
boson pairs {VV , V=W, Z) where one boson decays to a dijet final state . The data correspond to 
3.5 fb~^ of integrated luminosity oipp collisions at ^/s = 1.96 TeV collected by the CDF II detector 
at the Fermilab Tevatron. We observe 1516 ± 239(stat) ± 144(syst) diboson candidate events and 
measure a cross section a{j)p^VV+X) of 18.0 ± 2.8(stat) ± 2.4(syst) ± l.l(lumi) pb, in agreement 
with the expectations of the standard model. 
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PACS numbers: 14.80.Bn, 14.70.-e, 12.15.-y 



The production of heavy gauge boson pairs {WW, WZ, 
or ZZ) in pp coUisions has been observed in the fully 
leptonic final states at the Fermilab Tevatron collider [TJ 
[5]. Diboson production has not yet been conclusively 
observed in pp collisions in decay channels involving 
hadrons [3]; however, evidence for diboson decays into 
an Ivqcf'^ final state (/ = e, /u, r; q = u, d, s, c, b) has been 
recently presented by D0 collaboration [4]. 

Measurements of diboson production cross sections 
provide tests of the self-interactions of the gauge bosons. 
Deviations from the standard model (SM) prediction 
for the production rates could indicate new physics [5j. 
Furthermore, given that diboson production is topo- 
logically similar to associated Higgs boson produc- 
tion, pp^VH+X (V^W, Z), the analysis techniques de- 
scribed in this Letter are important for Higgs boson 
searches. 

Here we present the first observation at a hadron col- 
lider of diboson production with one boson decaying into 
leptons and the other into hadrons. The analysis is per- 
formed on a sample of events with large transverse mo- 
mentum imbalance (-^t) [ElIZ] and two jets whose invari- 
ant mass can be reconstructed. This signature is sensitive 
to not only Ivqq^'^ , but also to vvqq^'^ decays because we 
do not explicitly require presence of identified charged 
leptons. The limited dijet mass resolution results in a 
significant overlap of the W^qq' and the Z^qq dijet 
mass peaks, and therefore the combination of the three 
diboson signals, WW, WZ, and ZZ, is considered. 

We analyze a dataset of pp collisions corresponding 
to an integrated luminosity of 3.5 fb^^ collected with 
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the CDF n detector at the center-of-mass energy of 
1.96 TcV. The CDF H detector is described in detail 
elsewhere [Sj. Here the components that are relevant to 
this search are briefiy discussed. Surrounding the beam 
pipe, there is a tracking system consisting of a silicon 
microstrip detector, a cylindrical drift chamber, and a 
solenoid that provides a 1.4T magnetic field along the 
beam axis. The central and plug calorimeters, which 
respectively cover the pseudorapidity regions of |77|<1.1 
and 1.1<|77|<3.6 [B], surround the tracking system with 
a projective tower geometry. The calorimeters are com- 
posed of inner electromagnetic and outer hadronic sec- 
tions that consist of lead-scintillator and iron-scintillator, 
respectively. In the central region, the calorimeter con- 
sists of 48 modules, segmented into towers of granu- 
larity A?7 X A(j> « 0.1 X 0.26. The energy resolution 
of the central electromagnetic calorimeter for electrons 
is (j{Et)/Et = 13.5%/v^i;T(GeV) ® 1.5% [9], while 
the energy resolution of the central hadron calorime- 
ter for charged pions that do not interact in the elec- 
tromagnetic section is a{ET)/ET = hQ% / ^ EriGeSf) ® 
3% [TU], where Et is the transverse energy [B]. The 
plug electromagnetic (PEM) and plug hadron (PHA) 
calorimeters are identically segmented into 480 towers 
whose sizes are At; x A</) ~ 0.1 x 0.13 for |r7|<1.8 
to A?7 X A(^ ~ 0.6 X 0.26 for the towers centered at 
|?7|=3.3. The corresponding PEM and PHA energy reso- 
lutions are a{ET)/ET = 14.4%/ y^iJ^^GeV) ® 0.7% and 
a{ET)/ET = 74%/yi^T(GeV) ® 4%, respectively [U]. 
The wall hadron calorimeter (WHA) covers the gap in ac- 
ceptance between the central and plug hadron calorime- 
ters, corresponding to 0.7<|77|<1.3, with segmentation 
similar to that of the central calorimeter. The energy 
resolution of WHA is a{ET)/ET = 75% / V^MCeV) ® 
4% [in] for charged pions that do not interact in the 
electromagnetic section. A system of Cherenkov coun- 
ters [12], located around the beam pipe and inside the 
plug calorimeters, is used to measure the number of in- 
elastic pp collisions per bunch crossing and thereby the 
luminosity. 

The diboson signal [WW, WZ, and ZZ) is simu- 
lated using the pythia v6.2 Monte Carlo generator [13]. 
The most significant backgrounds to the diboson sig- 
nal are W{lv)+]eis, Z{vv)+]ets, and QCD multijet pro- 
duction (in the following referred to as multijet back- 
ground or MJB). Other less significant backgrounds in- 
clude Z(Z/)-f jets, tt, and single t-quark production. The 
W^-|-jets backgrounds are simulated using the fixed-order 
matrix element generator ALPGEN v2.1 [2] which is in- 
terfaced with PYTHIA v6.3 to simulate parton showering 
and fragmentation, the underlying event, and additional 
pp interactions in the same bunch crossing. The Z-|-jets 
and t-quark production processes are simulated with 
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PYTHIA v6.2. The detector response in all Monte Carlo 
samples is modeled by a GEANT-based CDF II detector 
simulation [15] . The MJB does not typically result in sig- 
natures of large intrinsic IpT- However, when jet energy is 
not measured accurately an event may be reconstructed 
with large Ifyr and pass the analysis selection criteria. Be- 
cause of the large multijet production rate, this can still 
be a significant background in a -^-l-jets based analysis. 
The MJB is determined from the data. All other back- 
ground predictions are normalized using next-to-leading- 
order (NLO) calculations for SM cross sections. 

The selection of signal events proceeds as follows: a 
set of ^T-based triggers select events with a variety 
of IfjT and jet requirements. All these triggers have 
benefited significantly from the calorimeter trigger up- 
grade completed in 2007 [16]. The majority (94%) of 
events satisfy the inclusive Iflr trigger which requires 
^T>45GeV. Jets are reconstructed in the calorimeter 
using the jetclu cone algorithm [17^ with a cone ra- 
dius of 0.4 in (?7, (j)) space. The energy measured by the 
calorimeter is corrected for effects that distort the true 
jet energy [18]. Such effects include the non- linear re- 
sponse of the calorimeter to particle energy, loss of en- 
ergy in uninstrumcntcd regions of the detector, specta- 
tor interactions, and energy radiated outside of the jet 
cone. We select events that have ^T>60GeV [H] and 
exactly two jets with £^T>25GeV and |77|<2.0. This 
ensures a trigger efficiency of 96%±2% on signal. We 
search for diboson production in the dijet mass range 
40 <Mjj<160 GeV/c^. The lower edge of the mass range 
is chosen to ensure that events are on the trigger effi- 
ciency plateau with respect to dijet mass. In addition 
to the requirements discussed above, the electromagnetic 
fraction of the total energy for each of the two jets is re- 
quired to be less than 90% to ensure that electrons and 
photons are not counted as jets. 

In order to suppress the MJB we use a Ifrp resolution 
model to distinguish true Ipr originating from undetected 
neutrinos from fake IpT due to jets that are not measured 
accurately. The I^t significance is a dimensionless quan- 
tity based on the energy resolution of the jets, on soft 
unclustered particles, and on the event topology; details 
of its definition can be found in Ref . [20 . The If^x signif- 
icance is typically low when ]pT arises from mismeasure- 
ment. In addition to having a small significance, the Iflj, 
will often be aligned with a jet. We select events with Ifrr 
significance larger than 4 and azimuthal angle between 
and nearest jet (A^^^) greater than 0.4 radians. 

Finally, we apply several requirements that suppress 
contamination due to cosmic-ray, beam-related, and 
other non-collision backgrounds. Events are required to 
have at least one reconstructed vertex formed by charged 
particle tracks. The transverse energies of all calorime- 
ter towers are calculated with respect to the z position 
of the primary vertex with the largest of associ- 



ated tracks. The electromagnetic fraction of the total 
event energy has to be larger than 30% in order to re- 
duce beam-related backgrounds. The arrival time of both 
leading jets as measured by the electromagnetic shower 
timing system [5T] has to be consistent with the pp col- 
lision time. The remaining non-collision background has 
a smooth Mjj distribution and accounts for less than 90 
of the 44,910 selected events. 

The shape and normalization of the MJB are deter- 
mined from the data. A vector, analogous to the 

calorimeter-based If^rp, is constructed from the vector 
sum of the transverse momenta of particles measured 
in the tracking system, and is largely uncorrelated to 
If^rp for events where jets are not reconstructed accu- 
rately. In the absence of IpT arising from mismeasure- 
mcnt in the calorimeter, the If^j, and ^j, will be aligned 
in most events. The MJB is expected to be the dominant 
background component at larger values of A0(_^y,^y). 
The dijet mass shape and normalization for the remain- 
ing MJB contribution in the sample is found by select- 
ing events with l^(f){Iprp,]/)rp)>\.Q and subtracting out the 
non-MJB backgrounds. The normalization is scaled up 
to account for the MJB contamination in the region 
A0(^j,, jij,)<1.0. The shape of the MJB is fit to an expo- 
nential in Mjj to derive a dijet mass template. The MJB 

shapes of Mjj and A(j){^j,,-^rp) distributions are verified 
with a large statistics MC sample. 

The signal extraction is performed using a minimiza- 
tion of the unbinned extended negative log likelihood 
with RooFiT program [55]. Three Mjj template distri- 
butions are used in the fit: the first is F-l-jets and t-quark 
production (in the following referred to as "electroweak" 
(EWK) backgrounds) and is taken from Monte Carlo sim- 
ulation; the second is the MJB template, where the slope 
and normalization are Gaussian constrained to their pre- 
viously measured values; the third template describes the 
signal. The signal shape is comprised of the WW, WZ, 
and ZZ distributions. This template is obtained from a 
Gaussian-f polynomial fit to the signal Monte Carlo sim- 
ulation where the mean and the width of the Gaussian 
distribution are linearly dependent on the jet energy scale 
(JES). 

To assess the effect of systematic uncertainties on 
the measurement, we address separately two classes of 
sources: those that affect the signal extraction proce- 
dure, and those that affect the signal acceptance in the 
cross section calculation. The signal extraction system- 
atic uncertainties come from uncertainties in signal and 
background shapes. The shape uncertainties take into 
account the effect of jet energy resolution (JER), JES, 
MJB shape, and the shape of the EWK background. The 
jet energy scale, and the shape and the normalization of 
MJB, are treated as nuisance parameters in the fit and 
Gaussian constrained to their independently measured 
values. These uncertainties are therefore accounted for 
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in the statistical uncertainty of the extraction. 

The shape uncertainty for the EWK background is 
determined by using 7+jets data [23] as an alternative 
background model in the Mjj fit. All major non-MJB 
backgrounds include a gauge boson accompanied by jets. 
There are similarities between the 7+jets and l/+jets 
production; however, due largely to the mass difference 
between the 7 and the W/Z, the kinematics is not iden- 
tical. To take this into account the 7+jets data are 
weighted by the ratio of the dijet mass distributions of the 
EWK background MC samples to 7+jets pythia sam- 
ple. We use these adjusted 7+jets data to determine a 
systematic uncertainty on the EWK Mjj template. Se- 
lection cuts applied to 7+jets events are not identical to 
those applied to the ^T+jets sample. For example, the 
Z decay into neutrinos will register as in the detector, 
while the photon Et will be measured in the calorimeter. 
For this reason we cut on the vector sum of the photon Et 
and any Ipr present in 7+jets events at 60 GeV, treating 
this sum as analogous to Ipr in V+jets events. A further 
consideration in the construction of the 7+jets template 
is the effect of events, as these events will cause 
a peak in the 7+jets dijet mass distribution. We sub- 
tract this contribution using the j+V pythia sample. 
Finally, we perform two signal extraction fits using the 
default EWK and 7+jets templates, respectively. The 
uncertainty due to the shape of the EWK background is 
then estimated as the difference in the results obtained 
from these two fits. The described method accounts for 
a combined effect of JES, JER, and modeling of jets in 
MC on the EWK Mjj template. 

The uncertainty associated with the JES is the domi- 
nant source of systematic uncertainty on the acceptance, 
and, therefore, the cross section. The JES affects sev- 
eral of the variables used in the event selection. The 
effect of the JES systematic uncertainty is quantified by 
varying the jet energies in the signal Monte Carlo simula- 
tion to account for the ±lcr variations of the JES. Other 
less significant sources of systematic uncertainty that af- 
fect the measured cross section are jet energy resolution, 
initial and final state radiation (ISR/FSR), and parton 
distribution functions (PDF). A summary of all sources 
of systematic uncertainty is presented in Table [T] 

The measured yields for signal and backgrounds are 
given in Table |lT] We extract 1516±239(stat) ±144(syst) 
signal events which correspond to a cross section of 
18.0 ± 2.8(stat) ± 2.4(syst) ± l.l(lumi) pb, in agreement 
with the SM prediction of 16.8+0.5 pb obtained using the 
MCFM v5.4 program [H] with CTEQ6.1M PDFs [25]. 
Based on the MC simulation, the acceptances for the 
WW, WZ, and ZZ production is 2.5%, 2.6%, and 2.9%, 
respectively. In the calculation of the combined diboson 
cross section, we assume that each signal process con- 
tributes proportionally to its predicted SM cross section: 
11.7 pb for WW, 3.6 pb for WZ, and 1.5 pb for ZZ. 

Figure [T] shows a comparison between the observed 



A(/)-^^ distribution and the MJB and EWK (sig- 
nal+background) components. This distribution pro- 
vides a strong consistency check on our MJB model. Fig- 
ure |2] shows the fit result and a comparison between the 
expected signal and data after background subtraction. 
We bin the data as in Fig. [2] and obtain a of 9.4 for 9 
degrees of freedom corresponding to a p- value of 40%. 

In summary, we use the ^+jets final state to measure 
the WW+WZ+ZZ cross section in pp collisions at -y/s — 
1.96 TeV to be 18.0±2.8(stat)±2.4(syst)±l.l(lumi)pb. 
This is consistent with the SM prediction of 16.8+0.5 pb. 
To assess the strength of the observed signal, the ef- 
fect of parameter variations due to all relevant sources 
of uncertainty are studied by comparing the likelihood 
of background-only fit with the full fit result, and con- 
verting the difference into significance numbers. We thus 
measure that the signal corresponds to a significance of at 
least 5.3 standard deviations from the background-only 
hypothesis. 

TABLE I: The systematic uncertainties and their effect on 
the number of extracted signal events, the acceptance, and 
the cross section. All systematics are added in quadrature. 





Systematic 


% uncert. 


Extraction 


EWK shape 
Resolution 


7.7 
5.6 




Total extraction 


9.5 




JES 


8.0 




JER 


0.7 


Acceptance 


]pT resolution model 1.0 
Trigger inefficiency 2.2 
ISR/FSR 2.5 
PDF 2.0 




Total acceptance 


9.0 




Luminosity 


5.9 




Total 


14.4 



TABLE II: Value of parameters in the model used to fit the 
Mjj distribution. 1516±239(stat) signal events are extracted 
from the 44,910 data events which pass our selection cuts. 
The jet energy scale is also extracted from the fit and agrees 
well with the default value (1.0) as measured from calibra- 
tions "18'. 

Parameter Fitted value 

Jet energy scale, JES 0.985 + 0.019 

Yield of EWK background events 36, 140 + 1230 
Yield of MJB events 7249 + 1130 

Yield of diboson candidates 1516 + 239 
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FIG. 1: Data compared with the sum of the predicted EWK 
and MJB backgrounds for the A(^^'|^ variable. The band rep- 
resents the total systematic uncertainty on the background. 
The measured signal is included here in the EWK contribu- 
tion. 
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FIG. 2: Top: Comparison between data and fitted back- 
ground only. The measured signal is shown unstacked. The 
band represents the systematic uncertainty due to the shape 
of EWK background as described in text. Bottom: Compar- 
ison of the diboson signal (solid line) with the background- 
subtracted data (points). The dashed lines represent the ±la 
statistical variations on the extracted signal. The gray band 
represents the systematic uncertainty due the EWK shape. 
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